Density of states of disordered graphene 
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We calculate the average single particle density of states in graphene with disorder due to impurity 
potentials. For unscreened short-ranged impurities, we use the non-self-consistent and self-consistent 
Born and T-matrix approximations to obtain the self-energy. Among these, only the self-consistent 
T- matrix approximation gives a non-zero density of states at the Dirac point. The density of states 
at the Dirac point is non-analytic in the impurity potential. For screened short-ranged and charged 
long-range impurity potentials, the density of states near the Dirac point typically increases in the 
presence of impurities, compared to that of the pure system. 
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I. INTRODUCTION 

The recent experimental realization of a single layer 
of carbon atoms arranged in a honey-comb lattice has 
prompted much excitement and activity in both the ex- 
perimental and theoretical physics communities^'^. Car- 
riers in graphene (both electrons and holes) have a lin- 
ear bare kinetic energy dispersion spectra around the 
K and K' points (the "Dirac points") of the Brillouin 
zone. The ability of experimentalists to tune the chemi- 
cal potential to lie above or below the Dirac point energy 
(by application of voltages to gates in close proximity to 
the graphene sheets) allows the carriers to be changed 
from electrons to holes in the same sample. This sets 
graphene apart from other two dimensional (2D) car- 
rier systems that have a parabolic dispersion relation, 
and typically have only one set of carriers, i.e. either 
electrons or holes. Another unique electronic property, 
the absence of back-scattering, has led to the specula- 
tion that carrier mobilities of 2D graphene monolayers 
(certainly at room temperature, but also at low temper- 
ature) could be made to be much higher than any other 
field-effect type device, suggesting great potential both 
for graphene to be the successor to Si-MOSFET (metal- 
oxide-semicondcutor field effect transistor) devices and 
for the discovery of new phenomena that normally ac- 
companies any significant increase in carrier mobility'^"^ . 
It is therefore of considerable fundamental and techno- 
logical interest to understand the electronic properties of 
graphene^. 

Graphene samples that are currently being fabricated 
are far from pure, based on the relatively low electronic 
mobilities compared to epitaxially grown modulation- 
doped two-dimensional electron gases (2DEGs) such as 
GaAs-AlGaAs quantum wells. It is therefore important 
to understand the effects of disorder on the properties of 
graphene. Disorder manifests itself in the finite lifetimes 
of electronic eigenstates of the pure system. In the pres- 
ence of scattering from an impurity potential that is not 
diagonal in the these eigenstates, the lifetime scattering 
rate of the eigenstate, 7, is non-zero and can be measured 



experimentally by fitting the line-shape of the low-field 
Shubnikov-de Haas (SdH) oscillations.'''^ The effect of 
disorder scattering on the SdH line shape is equivalent 
to increasing the sample temperature and one can there- 
fore measure this Dingle temperature {Tb) and relate it 
to the single-particle lifetime through /17 — iirkBTrt- To 
avoid potential confusion, we mention that the lifetime 
damping rate 7 discussed in this paper is not equal to 
the transport scattering rate which governs the electrical 
conductivity. The lifetime damping rate is the measure 
of the rate at which particles scatters out of an eigen- 
state, whereas the transport scattering rate is a mea- 
sure of the rate of current decay due to scattering out 
of an eigenstate. In normal 2DEGs, the transport scat- 
tering time can be much larger than the impurity in- 
duced lifetime h/2"f, particularly in high mobility mod- 
ulation doped 2D systems where the charged impurities 
are placed very far from the 2DEG®'^. Recently, the issue 
of transport scattering time versus impurity scattering 
lifetime in graphene has been discussed. 

The single-particle level broadening due to the impu- 
rity potential changes many of the physical properties of 
the system including the electronic density of states. ^^'^^ 
The electronic density of states is an important property 
which directly affects many experimentally measurable 
quantities such as the electrical conductivity, thermoelec- 
tric effects, and differential conductivity in tunneling ex- 
periments between graphene and scanning-tunneling mi- 
croscope tips or other electron gases. Changes in the 
density of states also modify the electron screening, ^"^ 
which is an important factor in the determination of var- 
ious properties of graphene. It is therefore imperative 
to take into account the effects of disorder on the den- 
sity of states, particularly since disorder is quite strong 
in currently available graphene samples. 

In the present work, we present calculations of the av- 
erage density of states of disordered graphene. This prob- 
lem has been investigated using various models and tech- 
niques, both analytical and numerical. ^'^'^ We take 
into account scattering effects from long-range and short 
range impurity potentials. We consider both unscreened 
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and screened short-ranged and screened charged impuri- 
ties, using the Born approximation. In addition, for un- 
screened short-ranged (USR) impurities, we go beyond 
the Born approximation and include self- consistent ef- 
fects. 

There is another class of disorder in graphene called 
"off-diagonal" or "random gauge potential" disorder, in 
which the hopping matrix elements of the electrons in the 
underlying honeycomb lattice are random. In this paper 
we do not consider in this type of disorder, which can 
result from height fluctuations (ripples) in the graphene 
sheet and lead to qualitatively different results from the 
ones presented in this paper. ^'^'-'^^'^^'^^ 

The rest of the paper is organized as follows: In Sec. II, 
we describe the approximation schemes that we use. 
Sees. Ill and IV deal with unscreened short range impuri- 
ties and screened short-range/charged impurities, respec- 
tively. In Sec. V, we compare our results to those from 
other workers, and we conclude in Sec. VI. 

II. APPROXIMATIONS FOR THE 
SELF-ENERGY 

The single-particle density of states for a translation- 
ally invariant 2DEG is given by^* 

^(^^) = -fE/(^MGA(fe,i;)] (1) 



For pure graphene systems, S(fc, E) = Q (ex- 
chiding electron electron and electron phonon inter- 
actions, which are not considered here), and hence 
Im[GA(A;, £^)] = —Tr6{E — Ek\)- Close to the Dirac points 
(which we choose to the the zero of energy), the disper- 
sion for graphene is (we use ft = 1 throughout this paper) 

Ek,x = Xvpk, (4) 

where A = -1-1 and —1 for the conduction and valence 
bands, respectively, k = \k\ is the wavevector with re- 
spect to the Dirac point, and Wi? is the Fermi velocity 
of graphene. Performing the fe-integration in Eq. (1) for 
the pure graphene case gives 

Do{E) = f\§9{\E\-E,), (5) 

where E,. = vpkc is the band energy cut-off. 

The average density of states for a disordered 2DEG 
can be obtained by averaging the Green's function over 



where g is the degeneracy factor (for graphene g = A due 
to valley and spin degeneracies) A is the band index, G 
is the retarded Green's function, and the fe-integration 
is over a single valley, which we assume to be a circle of 
radius \k\ = kc- G expressed in terms of the retarded 
self-energy T,\{k,ULi) is 



Gx{k, E) = [E- Ek,x - SA(fe, E) + iri\-\ (2) 



where Ek \ is the bare band energy of the state \k\) 
and rj is an infinitesimally small positive number. (In 
this paper, the Green's functions and self-energies are all 
assumed to be retarded.) Eqs. (1) and (2) show that if 
Ek^x and Y,x{k,E) are known, the density of states can 
be obtained in principle from 



(3) 

I 

impurity configurations. The averaging procedure gives 
a non-zero S which, in general, cannot be evaluated ex- 
actly. Various approximation schemes for S have there- 
fore been developed, four of which are described below. 

Bom Approximation — In the Born approximation, 
the self-energy is given by the Feynman diagram shown 
in Fig. l(a)^'', and the expression for the self-energy is 

/dk' 
(2^ 

xY,Go,x'{k',E)Fxx'{k,k') (6) 

A' 

where m is the impurity density, U{q) is the Fourier 
transform of the impurity potential, Go is the bare 
Green's function and Fxx' {k, k') is square of the the over- 
lap function between the part of the wavefunctions of 
|fcA) and |fc'A') that are periodic with the lattice (here 
A, A' are band indices). For graphene states near the 



dk -lTa[T.x{k,E)]+ri 
WF {E - Ek,x - RepA(fe, E)])^ + {lm[Ex{k, E)] - nf 
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(a) 



X 



matrix clement joining intcrvalley states is constant, in- 
clusion of intcrvalley scattering in our calculations is not 
difBcult.) This simplification allows us to obtain some 
analytic expressions for the self-energies in the approxi- 
mation schemes mentioned above. 




A. Self-energy for USR disorder 

For USR disorder, the four approximation schemes we 
use give self-energies that are independent of A and k. 



FIG. 1: Feynman diagrams for the (a) Born and (b) T-matrix 
approximations for the self-energy. The "x" , dotted line and 
line with arrow signify the impurity, impurity potential, and 
Green's function, respectively. The Green's functions are ei- 
ther bare or self-consistent. 



Dirac point, 



Fxy{k,k') = -{l + X\'cosl 



(7) 



where 6kk' is the angle between k and k' . 

Self-Consistent Born Approximation — The Feynman 
diagram for this self-energy is the same as Fig. 1(a), 
except that the bare Green's function is replaced by the 
full one. Consequently, the expression for Esb is the same 
as in Eq. (6), except with Gq replaced by G. 

T-matrix Approximation — The T-matrix approxima- 
tion is equivalent to the summation of Feynman diagrams 
shown Fig. 1(b). The expression for St is the same as 
in Eq. (6), except with U replaced by T, the T-matrix 
for an individual impurity. 

Self-consistent T-matrix Approximation — The Feyn- 
man diagram for this approximation is the same as in the 
T-matrix approximation, except that the bare Green's 
functions are replaced by full ones. The expression for 
Est is the same as in Eq. (6), except with U replaced by 
T, and Go replaced by G. 

When the potentials for the impurities are not all iden- 
tical (for example, in the case where there is a distribu- 
tion of distances of charged impurities from the graphene 
sheet), one averages \U{k — fe')|^ or |T(fe — fe')|^ over the 
impurities. 



III. UNSCREENED SHORT-RANGED (USR) 
DISORDER 



In the present context, short-ranged impurities are im- 
purities which result from localized structural defects in 

the honeycomb lattice, which are roughly on the length 
scale of the lattice constant. In this case, it is accept- 
able to approximate U{q) = Uo, a real constant, for in- 
travalley scattering processes. In this paper, we ignore 
the intcrvalley processes. (However, we note that if the 



1. Bom approximation 

The self-energy for graphene with USR scatterers in 
the Born approximation, using U{q) = Uq, Eqs. (4) and 
(7) in Eq. (6), is 
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As a function of complex (, i?o(C) is real and positive 

(negative) along the real axis from Ec to oo {—Ec to — oo). 
Furthermore, it has a branch cut in on the real axis of ^ in 
between —Ec and Ec, so that for —Ec<E (real) < Ec, 



Ho{E±ir]) = -(27r)^^ i^Eln 

Fig. 2 shows Hq [E + irj) . 
For real E and \E\ < E^ 
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(10) 



The Born approximation damping rate for state |feA) is 
jB{k) = -2lm[i:''^"'\Ek,x)] = iBVpk. 



2. Self- consistent Bom approximation 

The self-energy Eg"^"^^ for unscreened short-ranged 
scatterers in the self-consistent Born approximation is 
given by the self-consistent equation 

^(usr) , 



S^B n^)=7B^SB(i? + ir/); 



(11a) 



i^SB(C) = 
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This shows that if < then T.'^sbH^) 



j^^sr)^^-j j'gjj.pgp^ possibly around E ~ ±Ec, which is 
usuaUy not experimentaUy relevant). 



3. T-matrix approximation 



FIG. 2; Real and imaginary parts of Ho(E + irf) — 

^b"\E)/^b [see Eq. (9)], where T.'g'''\E) is the self-energy 
for USR impurities in the Born approximation. 



In general the impurity averaged T-matrix for a poten- 
tial U{q) is 



Tko\o,k\{E) -rii^ 



n=l 



n H /-^C/(fc,_i-fcOGo,A.(fe.,S)FA._,,A,(fc^-i,fcO) C^(fe^-fc)^A„,A(A;„,fc) . 

02) 



In this approximation, the self energy is 
j:TAk,E)^Tkx,kx{E)- 



As in the self-consistent Born approximation, this shows 
that if |47''(-S)I < then 1:''st\e) « 

(13) 

(except possibly around E — ±Ec)- 



If U{q) = Uq, a constant, the term in the square paren- 
theses in Eq. (12) is Uq^'^HS{E + ir]) / {2vj,)'' . The sum 
then gives 



^BHo{E + if]) 



4 



(14) 



4- Self- consistent T-matrix approximation 

As in self-consistent Born approximation, the Hq{E) in 
Eq. (14) is replaced by the self-consistent Hst{E), giving 



-^(usr) 



{E) = 



^bHst{E + irj) 



1-^Hst{E + i7^) 



(15a) 
(15b) 



B. Density of states for fe-independent S 



When the self-energy is A;- and A- independent, the den- 
sity of states for graphene can be calculated analytically 
from Eq. (3). The result is 



D{E) 




tan 



Er-n 



tan 



Er + n 



2 tan 



(16) 
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where T{E) = -Im[I](^)] and 0(E) =E- Re[S(£;)]. 



gives 



C. Density of states at the Dirac point 

There has been considerable interest in the minimum 
DC electrical conductivity of disordered graphene as the 
Fermi energy moves through the Dirac point. ^^'^^ There 
is still no consensus on whether the minimum conductiv- 
ity is a universal value or not. Since the electrical con- 
ductivity is directly proportional to the density of states 
at the Fermi energy, it is important to be able to deter- 
mine the density of states at the Dirac point of disordered 
graphene. Because the minimum conductivity is non-zero 
as the Fermi energy passes through the Dirac point, the 
density of states should be nonzero. 

Eq. (3) shows that if T.{E) ^ when E ^ 0, then 
the density of states at the Dirac point D{0) = 0. (Note 
that we have not included the term that is first order in 
the impurity potential in our self-energy. Since this first- 
order term merely rigidly shifts the band by an amount 
riiUo, ignoring this term is equivalent to shifting the zero 
of the energy by —riiUo, and hence the Dirac point is still 
at = 0.) A non-zero density of states at the Dirac point 
depends on a non-zero Im[Yj{E = 0)]. Since Hq{E — > 
0) = 0, it is clear from Eqs. (8a) and (14) that the Born 
and T-matrix approximations give zero density of states 
at the Dirac point. 

For the self-consistent Born approximation, Eq. (lib) 
can be rewritten as 



-i(usr) 



(E) 



1 + 



2-K 



7b In 1 - 



(17) 



r^(o) ^ Im 



Er 



r^(o) 



= 27Br^(0)ln 



Er 



r«,(o). 



which implies that 



,(0) =£;eexp( - — 



(20) 



(21) 



Note that the result is non- analytic in Ua. Inserting 
Eq. (21) into the imaginary part of the denominator of 
Eq. (19) (which we had assumed to be much smaller 
than 1 in magnitude) gives the self-consistent criterion 
exp(-l/27B) < riiUo/Ec for the validity of Eq. (21). 
Substituting this into Eq. (16) gives an average density 
of states at the Dirac point for weak scattering of approx- 
imately 



27i 



(22) 



15,19,27 



Similar results to Eq. (22) have been reported 
studies of disordered systems of fermions with linear dis- 
persions using other methods. 

We mention that our calculation of the graphene den- 
sity of states at the Dirac point should only be consid- 
ered as demonstrative since electron-electron interaction 
effects are crucial^^ at the Dirac point, and the undoped 
graphene system is not a simple Fermi liquid at the Dirac 
point. 



IV. SCREENED SHORT-RANGED AND 
CHARGED IMPURITIES 



which shows that Sg^'"_^(E ^ 0) = 0, and therefore the 
self-consistent Born approximation also gives -D(O) = 0. 

In the case of the self-consistent T-matrix approxima- 
tion, re-writing Eq (15a) as 



E^7^(E)+niC/o = 



UiUo 



l-Ho{E-^^^^'\E))Uo/{2vj- 



(18) 

setting E = 0, using Eq. (8c) and taking the imaginary 
parts of both sides of this equation gives 



Im[E(0)] = r(0) = Im 



n,UQ 



• WO) 



Ei 



(19) 



In the weak scattering limit, when 73 <C 1, the imag- 
inary term in the denominator of Eq. (19) is much less 
than one (we check for self-consistency later), and this 



Free carriers will move to screen a bare impurity po- 
tential Vciiq), resulting in a screened interaction U{q) = 
Vci{q)/e{q), where e{q) is the static dielectric function. 
The e{q) results in an g-dependent effective electron- 
impurity potential U , even in the case of short-ranged 
(q-indepcndent) bare impurity potentials. This makes 
the calculations much more involved than in the USR 
case. Therefore, in this paper, we limit our investiga- 
tion of 5-dcpendcnt screened potentials to the level of 
the Born approximation. 

For the dielectric function e{q), we use the random 
phase approximation (RPA) for dielectric function appro- 
priate for graphene, given by e{q) — l — Vc{q)IiQ{q) where 
Vc{q) = /{nq) is the two-dimensional Fourier trans- 
form of the Coulomb potential (k is the dielectric con- 
stant of the surrounding material), and no(<7) is the static 
irreducible RPA polarizability for graphene. We use 
y&iil) = 14(?) for charged impurities, and 14i(9) =Uq, & 
constant, for short-range point defect scatterers. 

We first look at the density of states at the Fermi sur- 
face; i.e., at energy E = XkpVF- To obtain this, we 
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calculate the single-particle lifetime damping rate 7 in 
the Born approximation, which is given by 



lx{k) 



-2Im[SA(fc,SfeA)] 



Hi k 
2t: vp 



{\Ve^{q)?) 



(1 + COS 61), (23) 



where q — 2ks\n{9/2). Then, assuming that 
lm\Tix{k,E)] is relatively constant for k close to kp, 
we substitute \"i\{kF) for —lva\^x{k,EF)] into Eq. (3), 
which gives Eq. (16) with T = ^Jxikp). 

We assume that the charged or neutral impurities are 
distributed completely at random on the surface of the 
insulating substrate on which the graphene layer lies, the 
areal density for the charged and neutral impurities is 
riic and riis, respectively, and the density of carriers in 
the graphene layer is n = kp/w, where fc^? is the Fermi 
wavevector relative to the Dirac point. (This relationship 
between n and kp takes into account the spin and valley 
degeneracy gs = 2 and g^, = 2.) We use the RPA screen- 
ing function at T = 0,^^ to obtain the effective impurity 
potential. The key dimensionless parameter that quan- 
tifies the screening strength is rg = e^/^nvp), which is 
corresponding to the interaction strength parameter of a 
normal 2D system (i.e. the ratio of potential energy to ki- 
netic energy). The Born approximation lifetime damping 
rates for screened charged-impurities 7c and (5-correlated 
neutral impurities 75 at /cf are 



isikp) = Is{2rg). 



In these equations, 



Ic{x) ^x- — +x^f{x); 



(24a) 
(24b) 

(25a) 



Isix) = 
where 



vr 
4 



3x(l-^^ +x{3x^ ^2)f{x), (25b) 



l+Vl-a^ 



-1 1 



for X < 1; 

for X — I] 
for X > 1. 



(26) 



In Fig. 3 we show the calculated damping rates scaled 
by \Ep\/ = kpvp as a function the interaction param- 
eter Ts. For <C 1 we have jc/Ep ft! nicrs/2n and 
'ys/Ep ~ "SbI^- For :s> 1 we have ^c/Ep ~ nnic/lGn 
and js/Ep « 7B/(27rrs). Thus, for small (large) the 
damping rate due to the short-ranged impurity domi- 
nates over that due to the long-ranged charged impu- 
rity. On the other hand, since ^c{kp) oc kp^ oc and 
"fsikp) (X kp (X ni [Eq. (24)], in the low (high) carrier 
density limit the lifetime damping of single particle states 
at the Fermi surface is dominated by charged impurity 




FIG. 3: Calculated damping rates scaled by Fermi energy 
^/Ef as a function of Vs- Graphene on a Si02 (air) sub- 
strate has an ~ 0.7 (2). Solid lines indicate damping 
rates (7c) due to charged impurities with an impurity den- 
'^cm"^ for different electron densities n = 1, 10, 
cm~^ (from top to bottom), respectively. Dashed 



sity riic = 10 
50x10 

line indicate the damping rate (75 
purity with impurity density nu = 

strength Uo ~1 KeV A'^, which correspond to 7s — 0.11 
Note 'ys/Ep is independent on the electron density. 



due to short-ranged im- 
10^^ cm"^ and potential 



(short-ranged impurity) scattering. The crossover takes 
place around a density 



n,c nvp Ic{2rs) 
^s^Is{2rs)' 



(27) 



Using r = "f\{kp)/2 in Eq. (16) gives (assuming Ep, 
r « E,) 



D{Ep) « Do{Ep) 



1 1 

- H — tan 

2 TT 



In 



\Ep 



2'k\Ep 

For T/\Ep\ < 1, this gives 
D{Ep)^Do{Ep)\l + -- ^ 



E'j. 



El + V^ 



(28) 



and for \Ep\/T < 1 
D{Ep) « D^{Ep) 



TT \Et. 



\Et. 



In 



E, 
\Ep 



(29) 



ttF 



9s9v__£_ 
27r2 vl 



(30) 

We can apply Eq. (29) for short-ranged impurity scat- 
tering and for charged impurity scattering in high carrier 
density limits, and Eq. (30) for charged impurity scat- 
tering in low density limits. Taking the limit \Ep\ — > 
in Eq. (30), it appears that for the case of screened 
Charged impurities, we obtain a finite density of states 
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CO/ Ep 

FIG. 4: (a) Imaginary and (b) real parts of self energy of a dis- 
ordered graphene ai k = kp for screened Coulomb scattering 
potential (solid lines) and for screened neutral short-ranged 
scattering potential (dashed lines). 



at the Dirac point with the Born approximation (since 
r oc 7c(fc_F) oc kp^). However, recall that in deriving 
Eqs. (29) and (30), we have assumed that T,{k,Ep) is 
constant with respect to fc, which is not necessarily the 
case at the Dirac point. 

In general, the damping rate (or the imaginary part of 
the self energy) is a function of energy and wave vector 
rather than a constant. From Eq. (6), we calculated the 
self-energy of disordered graphene. Fig. 4 we show the 
self-energy of a conduction band electron (A = -1-1) for 
both screened Coulomb scattering potential and screened 
neutral short-ranged potential. For Coulomb scatter- 
ers we use the impurity density riic = 10^^ cm~^, and 
for neutral short-ranged scatterers the impurity density 
riis = 10^^ cm~^ and potential strength Ua =1 KeV 
A^. The self-energies in the valence band (A — — 1) 
are related to the self-energy in the conduction band 
by ReS+(k,w) = -ReE_(k, -w) and ImS+(k,w) = 
ImE_(k, — w). As 0, ~lmS x{kF,uj) —>■ \uj\, and 

Iie'E,\{kp,LL!) a;ln|aj| for both scattering potentials. 
However, for large value of the asymptotic behaviors 
are different, that is, as \uj\ — > oo —lmT,{kp,Lu) oc |cl;|^"'^ 
for Coulomb scattering potential and — ImI](A:i?, cj) oc \uj\ 
for short-ranged potential. Note that by using only 
the (non-self-consistent) Born approximation in this sec- 
tion, we assume weak scattering and ignoring multiple- 
scattering events in calculated S(fc,w). Therefore, the 




co/Ep 

FIG. 5: The density of states in the presence of impurity (a) 
for screened Coulomb potential and (b) for screened short- 
ranged potential. In (a) we use the charged impurity densities 
TLic = 0, 1, 5 X 10^^ cm~^ (from bottom to top), and in (b) 
the short-ranged impurity density = 0, 0.5, 1 x 10^'^ cm~'^ 
(from bottom to top) and potential strength Uo =1 KeV A^. 



results are unreliable in the strong disorder limit (i.e., 
when Ti(k,uj) is modified significantly from its lowest- 
order form). 

In Fig. 5 the density of states in the presence of im- 
purity is shown for different impurity densities. In Fig. 
5(a) we show the density of states for Coulomb impurity 
potential with impurity densities, riic — 0, 10^^ cm~^, 
5 X 10^^ cm~^ (from bottom to top), and in Fig. 5(b) we 
show the density of states for neutral short-ranged impu- 
rity potential with densities, riid = 0, 5 x 10^^ cm~^, 10^^ 
cm^^ (from bottom to top) and potential strength Uq —1 
KeV A^ . The calculated density of states is normalized 
by Dq{Ef) — (gsf/t,/27r)£^F/7^. The density of states is 
enhanced near Dirac point {E — 0), but as \E\ ^ it 
goes zero as D{E) jE'l In for both types of screened 
impurity scattering. Based on the results of Section HI, 
we expect that this result is an artifact of the Born ap- 
proximation, and that the density of states should in fact 
be non-zero. The enhancement of the density of states 
can be explained as follows. In normal 2D system with 
finite disorder the band edge i?odgo,o of the pure conduc- 
tion (valence) band is shifted to i?cdgc.imp < (> 0) and 
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a band tail forms below (above) the band edge of a pure 
system. Thus, the density of states in the presence of 
impurities is reduced for £" > (i? < 0) because the 
states have been shifted by the impurity potential into 
the band tail. However, for graphene since the conduc- 
tion band and the valence band meet at the Dirac point 
the band tail (or shift of band edge) cannot be formed, 
which gives rise to enhancement of density of states near 
Dirac point. 

Before concluding we point out that our perturbative 
calculation of the graphene density of states assumes that 
the system remains homogeneous in the presence of impu- 
rities. It is, however, belicved^'^^ that graphene carriers 
develop strong density inhomogeneous (i.e. electron-hole 
puddle) at low enough carrier densities in the presence of 
charged impurities due to the breakdown of linear scat- 
tering. In such an inhomogeneous low-density regime 
close to the Dirac point, our homogeneous perturbative 
calculation does not apply. 

V. COMPARISONS TO OTHER WORKS 

In this section, we compare and contrast our model of 
disorder and results to other works in the field. 

Peres et al.^^ studied the effect of disorder in graphene 
by considering the effect of vacancies on the honeycomb 
lattice. For a finite density of vacancies, they found that 
the density of states at the Dirac point is zero for the 
"full Born approximation" (equivalent to our T-matrix 
approximation) and non-zero for the "full self-consistent 
Born approximation" (equivalent to our self-consistent T- 
matrix approximation). Our results are consistent with 
theirs, even though the regimes that are studied are dif- 
ferent. Vacancies correspond to the limit where the im- 
purity potential ?7o — > oo, whereas this work is more 
concerned with the weak impurity-scattering limit. 

Pereira at al}'^ considered, among several different 
models of disorder, both vacancies and randomness in the 
on-site energy of the honeycomb lattice. They numeri- 
cally calculated the density of states for these models of 
disorder. For compensated vacancies (same density of va- 
cancies in both sub-lattices of the honeycomb structure) 
they found that the density of states increased around 
the Dirac point. (Rcf. 17 also studied the case of un- 
compensated vacancies, but that has no analogue in our 
model of disorder.) For the case of random on-site impu- 
rity potential, they find that "there is a marked increase 
in the DOS (density of states) at Ed (the Dirac point)" 
and "the DOS becomes finite at Eo with increasing con- 
centration" of impurities. Our self-consistent T-matrix 
approximation result is consistent with their numerical 
results, although it should be mentioned again that it 
strictly does not apply to the case of vacancies. 



Wu et al. numerically investigated the average den- 
sity of states of graphene for the case of on-site disorder. 
They found that for weak disorder, the density of states 
at the Dirac point increased with both increasing den- 
sity of impurities and strength of the disorder. (In their 
work, they did not absorb the shift in the band due to 
the impurity potential in their definition of the energy, 
so the Dirac point had a shift Ed = xv where x is the 
concentration of impurities and v is the on-site impu- 
rity energy.) Their numerical results for the minimum in 
the average density of states (Fig. 3(b) in Ref. 20) seem 
to indicate a non-linear dependence of the value of the 
minimum as a function of the strength of the disorder 
potential, and is at least not inconsistent with Eq. (22). 

The issue of the effect of screening of the impurity in- 
teractions on the density of states discussed in Section IV, 
to the best of our knowledge has not yet been treated in 
the literature. Qualitatively, the effect of screened im- 
purities away from the Dirac point is the increase the 
density of states, and is consistent with numerical results 
for random on-site disorder. (Since our treatment of 
screened impurities is at the level of the Born approx- 
imation, we do not obtain either a non-zero density of 
states at the Dirac point nor resonances in the density of 
states.17'18.20) 



VI. CONCLUSION 

We have calculated the density of states for disordered 
graphene. In the case of unscreened short-ranged impuri- 
ties, we utilized the non-self-consistcnt and self-consistent 
Born and T-matrix approximations to calculate the self- 
energy. Among these, only the self-consistent T-matrix 
approximation gave a non-zero density of states at the 
Dirac point, and the density of states is a non-analytic 
function of the impurity potential. We investigated the 
density of states in the case of screened short-ranged and 
charged impurity potentials at the level of the Born ap- 
proximation. We find that, unlike the case of parabolic 
band 2DEGs, in graphene near the band-edge {i.e., the 
Dirac point) the density of states is enhanced by impu- 
rities instead of being suppressed. At very low carrier 
densities, however, graphene develops strong carrier den- 
sity inhomogeneity in the presence of charged impurities, 
an effect not captured by the homogeneous many-body 
theory in our description. 
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